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ABSTRACT 

The prompt optical flash recently detected accompanying GRB990123 suggests 
that, for at least some GRBs, 7-ray emission is accompanied by prompt optical-UV 
emission with luminosity L{1 — 7.5 eV) ~ 1 x 10^^(Ar2/47r)erg/s, where AQ is the solid 
angle into which 7-ray and optical-UV emission is beamed. Such an optical-UV flash 
can destroy dust in the beam by sublimation out to an appreciable distance, ~ 10 pc, 
and may clear the dust out of as much as ~ lO'^(Ail/47r)M0 of molecular cloud 
material on an apparent time scale of ~ ten seconds. Detection of time dependent 
extinction on this time scale would therefore provide strong constraints on the GRB 
source environment. Dust destruction implies that existing, or future, observations of 
not-heavily-reddened fireballs are not inconsistent with GRBs being associated with 
star forming regions. In this case, however, if 7-ray emission is highly beamed, the 
expanding fireball would become reddened on a ~1 week time scale. 

If the optical depth due to dust beyond « 8 pc from the GRB is 0.2 ^ ^ 2, most 
of the UV flash energy is converted to infra-red, A ~ l//m, radiation with luminosity 
LiR ^ lO^^erg/s extending over an apparent duration of « 20(1 + z)(AQ/0.01) day. 
Dust infra-red emission may already have been observed in GRB970228 and 
GRB980326, and may possibly explain their unusual late time behavior. 



Subject headings: Gamma Rays: Bursts- ISM: Dust, Extinction 



Introduction 



A prompt, 9-th magnitude, optical flash has recently been detected ( |Akerlof et al. 1999| ) 
accompanying the gamma ray burst (GRB) GRB990123. The most natural explanation of this 
flash is emission from a reverse shock propagating into fireball ejecta shortly after it interacts with 
surrounding gas ( |Sari &: Piran 1999b , Meszaros Sz Rees 1999 ). Although optical-UV emission 
from the reverse shock accompanying, or following shortly after, the 7-ray emission has been 
predicted (Meszaros, Rees & Papathanassiou 1994, Meszaros & Rees 1997, Panaitescu Sz Meszaros] 



T998i gari fc Piran 1999a| ) based on the simplest fireball models for gamma-ray bursts (GRBs), 
the intensity of optical emission could not have been reliably predicted due to the uncertainty in 
reverse shock parameters. The observations of GRB990123 suggest that the electron and magnetic 
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field energy fractions in the reverse sliock are similar to those in the forward shock propagating 
into the surrounding gas and producing the long term afterglow. This in turn implies that strong 
optical flashes accompanying 7-ray emission is a generic GRB characteristic. This is consistent 
with previous non-detection of optical flashes ( Park et al. 1997| ) given the wide GRB luminosity 



function ( [Krumholtz, Thorsett fc Harrison T998| , |Mau k MoT998| , [Hogg k Fruchter 19981 ) and the 



fact that GRB990123 is in the top 0.3% of the BATSE brightness distribution ( [Kouveliotou et 'al 



19991) 



There is evidence in several GRB afterglows for significant dust extinction, which may 



imply an association of GRBs with star forming regions ( [Kulkarni et al. 1998| , [Metzger et al 



1997 , Paczyiiski 1998 ). It has been shown that if GRBs indeed reside in such environment, the 
ionizing X-ray and UV afterglow radiation may lead to time dependent (on hour time scale) 
absorption ( perna &: Loeb 1998] ) and emission ( IGhisellini et al. 1998| , [Bottcher et al. 1998| ) line 



features. If H2 is present near the GRB, conspicuous 1110-1650Aabsorption will be imprinted 
on the spectrum, followed by UV fluorescence ( praine 1999b| ). On longer time scales, up to 
~ 10^ yr, GRB photo-ionization may lead to indicative recombination line features, which may 
allow identification of GRB remnants in nearby galaxies ( |Loeb &: Perna 1998 , Perna, Raymond 



fc Loeb 1999| ). Here, we discuss dust sublimation by the optical- UV flash accompanying the 



GRB. Since GRB990123 was an exceptionally bright burst, with exceptionally high intrinsic 7-ray 
luminosity, we first discuss in §2 the optical-UV flash itself to obtain the scaling of optical-UV 
luminosity with burst energy, and to estimate the prompt optical-UV emission at energies above 
the 2(1 + z) eV energy of the prompt emission observed by Akerlof et al. We show that if electron 
and magnetic field energy fractions in the reverse shock are similar among different bursts, then 
an optical-UV flash of luminosity Luv ~ 3 x 10^^(Ar2/47r)erg/s, where Ail is the beaming solid 
angle, is expected for a typical GRB. Dust destruction by thermal sublimation is discussed in §3, 
and other possible destruction mechanism are considered in §4. It should be emphasized that the 
analysis of the physics of dust destruction is independent of the model for optical-UV emission, 
and therefore of the discussion presented in §2. However, we present numerical results for dust 
destruction based on the characteristic optical-UV luminosity derived in §2. In §5 we discuss dust 
infra-red emission. The implications of our results are discussed in §6. 



2. The Prompt Optical-UV Luminosity 
2.1. Fireball Dynamics 



In fireball models of GRBs (see Piran 199£ for a recent review), the energy released by an 



explosion is converted to kinetic energy of a thin baryonic shell expanding at an ultra-relativistic 
speed. After producing the GRB, the shell impacts on surrounding gas, driving an ultra-relativistic 
shock into the ambient medium. After a short transition phase, the expanding blast wave 
approaches a self-similar behavior ( Blandford fc McKee 197^ ), where the expansion Lorentz factor 
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decreases with radius as F oc r~'^l'^ . The initial interaction of fireball ejecta with surrounding gas 
produces a reverse shock which propagates into and decelerates the fireball ejecta. Transition to 
self-similar behavior occurs on a time scale comparable to the reverse shock crossing time of the 
ejecta. 

The long term afterglow is produced by the forward, expanding shock that propagates 
into the surrounding gas. This shock continuously heats fresh gas and accelerates relativistic 
electrons, which produce the observed radiation through synchrotron emission. The most natural 
explanation of the optical flash is that it is due to synchrotron emission of electrons which are 
accelerated by the reverse shock, during the transition to self-similar behavior. Once the reverse 
shock crosses the ejecta, the ejecta expand and cool adiabatically. Thus, emission from the fireball 
ejecta is suppressed after the transition to self-similar expansion. 

Since the optical fiash is produced when the reverse shock crosses the ejecta, the plasma 
emitting the radiation expands with a Lorentz factor which is close to that given by the 
Blandford-McKee self-similar solution, F = (17-E'/167rnmpC^)^/^r~'^/^, where E is the fireball 
energy and n is the surrounding gas number density. The characteristic time over which radiation 
emitted by the fireball at radius r is observed by a distant observer is At ~ r/4r^c ( |Waxman 
1997cj ). The plasma Lorentz factor during optical fiash emission 



\ 10247rnmpC^ I \l + zj \ no J \1 + z 

where E = W^^E^^ erg, At°^- = lOAtf"- s is the observed duration. At"""- = (1 + z)At, and 
n = Iuq cm~^. 

Transition to self-similar expansion occurs on time scale At°^' comparable to the longer of 
the two time scales set by the initial conditions: the (observed) GRB duration AtQRB and the 
(observed) time Atr at which the self-similar Lorentz factor equals the original ejecta Lorentz 
factor Ti, r(At°'^- = Atr) = Fi. That is, 



At"^- = max 



1 + z fE^sV^^ fIi.Y^^^ 



2.5 V no / \300 / 



(2) 



During the transition, the unshocked fireball ejecta propagate at the original expansion Lorentz 
factor, Fj > F, and the Lorentz factor of plasma shocked by the reverse shock in the rest frame of 
the unshocked ejecta is ~ Fj/F. If At"^' ~ AtcRB ^ Atr then Fj/F 1, the reverse shock is 
relativistic, and the Lorentz factor associated with the random motion of protons in the reverse 
shock is ~ Fi/F. If At°^' ~ Atr ^ ^^grb then Fj/F ~ 1, and the reverse shock is not 
relativistic. Nevertheless, the following argument suggests that the reverse shock speed is not far 
below c, and that the protons are therefore heated to relativistic energy, 7^ — 1 ~ 1. The comoving 
time, measured in the fireball ejecta frame prior to deceleration, is tco. — ^/FjC. The expansion 
Lorentz factor is expected to vary across the ejecta, AFj/Fj ~ 1. Such variation would lead to 
expansion of the ejecta, in the comoving frame, at relativistic speed. Thus, at the deceleration 
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radius, tco. — TjAt, the ejecta width exceeds ~ ctco. — TjcAt. Since the reverse shock should cross 
the ejecta over a deceleration time scale, ~ TAt, the reverse shock speed must be close to c. We 
therefore conclude that the Lorentz factor associated with the random motion of protons in the 
reverse shock is approximately given by 7^ - 1 ~ Ti/r for both Ti/T ~ 1 and Fi/F > 1. For 
protons shocked by the forward shock 7^ ~ F, and therefore the ratio between thermal energy per 
proton in the reverse and forward shocks is (7^ — l)/(7^ — 1) ~ Tj/r^. Below we use this relation 
to derive the emission characteristics of the reverse shock by scaling the exact analytic results 



given for the forward shock emission in Gruzinov fc Waxman (1999) 



2.2. Optical-UV Emission 

If the fraction of thermal energy carried by electrons, ^e, and magnetic field, ^5, is similar 
in the forward and reverse shocks, then the frequency of peak synchrotron emission from the 
reverse shock is smaller than that of the forward shock by a factor ~ (Fj/r^)^. This is due 
to the fact that the energy density behind the reverse and forward shocks are similar, so that 
similar implies similar magnetic field strength in both regions, while similar implies 
le/le - (7^ - l)/(7^ - 1) ~ Fi/F2. Using Eq. (10) of [Gruzinov fc Waxman (1999)| for the 



forward shock peak frequency, we find that the reverse shock emission peaks at a frequency 



^ 1.7 X 10^^ {^Yf i^V^' ( I^Vnf Hz (3) 

^ 0.2J vo.oiy v3ooy ° ^ ^ 



(measured at the GRB redshift). Here and are the values relevant for the reverse shock, 
which in general may differ from those of the forward shock. The numerical values we have used 
are those characteristic of the forward shock (Waxman 1997a, Waxman 1997b, Wijers & Galama 



1998i [Granot et al. 1998|) 



As demonstrated below, the cooling time of electrons in the reverse shock, radiating in 
the optical-UV range, is long compared to the fireball expansion time. In this case, the peak 
synchrotron intensity is proportional to the product of magnetic field strength and number of 
radiating electrons (and independent of the electron Lorentz factor). The number of radiating 
electrons in the reverse shock is larger than in the forward shock by a factor F^/Fj. This can be 
deduced from the following considerations. The proton random Lorentz factor in the reverse shock 
is F^/Fj times smaller than that in the forward shock. Since the energy density in both regions is 
similar, the density of protons, and therefore of electrons, in the reverse shock is higher than that 
in the forward shock by the same factor. In addition, the width of the shocked fireball ejecta is 
similar to that of the shell of shocked surrounding gas, since both shocks propagate relativistically 
in the shocked plasma frame. Thus, if is similar in the reverse and forward shock, the peak 
synchrotron intensity fm is higher in the reverse shock by a factor F^/Fj. Using Eq. (11) of 
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Gruzinov fc Waxman (1999) we find for the reverse shock 



for a flat universe with zero cosmological constant, and Hq = 65/i65km/s Mpc. 

The decay of GRB990123 optical flux implies an electron energy distribution dNe/dE^ oc E~p 
with p ~ 2 ( Meszaros fc Rees 19991 ), for which the intensity at z/ > is oc i/"-*^/^. Thus, the 



optical [z^ = (5(1 + z) X 10 Hz] intensity is 



For the parameters of GRB990123, At°^' ~ 60 s and E ~ 10^^ erg (based on the 7-ray fluence 
~ 3 X 10~^erg cm~^), we obtain /y « 2 Jy, approximately twice the peak observed flux. Thus, 
the observed optical flash can be naturally explained by the simplest fireball model, provided the 
reverse shock parameters and are similar to those implied for the forward shock by afterglow 
observations, ~ 0.2 and ~ 0.01. If this is typical, than for a typical GRB, with E ~ 10^^ erg 
and At°^- « 10 s, we find 

where is the frequency at the redshift of the GRB. The form (|6|) is valid for frequencies u larger 
than the peak frequency z^m, and smaller than the frequency fc above which emission is dominated 
by electrons for which the cooling time is shorter than the dynamic time. Since the energy density 
in the reverse and forward shock regions is similar, fc (at the GRB redshift) is given by 

At frequencies u > Uc the spectrum steepens to L^, oc z^~^. 

Strong optical-UV emission requires, similar to GRB 7-ray production, large initial Lorentz 
factor, Fj > 100: Eq. (^) implies that At oc F^ for Fj ^ 300, and the optical-UV luminosity 
given by Eq. (^) drops as L,^ T?. This implies that the plasma emitting the optical-UV photons 
must also be expanding at high Lorentz factor, F > 100 [see Eqs. (§,|ll)], and strong optical-UV 
emission may be confined, similar to 7-rays, to a cone around the line of sight of opening angle 
1 ^ ^7 > 1/F ~ 0.01. This may be the case if, e.g., the fireball is a jet of opening angle 9j = 9^. 
A jet of finite opening angle expands as if it were a conical section of a spherical fireball, as long 
as Oj > 1/F. Thus, the analysis presented above is valid for a jet-like fireball. In this case, the 
energy E in the above equations should be understood as the energy that the fireball would have 
carried if it were spherically symmetric, and the optical-UV, as well as 7-ray, emission is confined 
to a small solid angle AQ = ttOj with optical-UV luminosity given by AilL^. 
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It should be noted that prompt optical emission has been observed for only 1 GRB to date, 
and therefore may not be typical of GRBs. Indeed, fireballs in low density environments with 
n <C lcm~^ would not be expected to produce strong prompt emission [see Eq. (^)]. However, we 
note from Eq. (^) that optical flashes are expected for typical GRBs with S53 1 and no 1, 
and should not be limited to unusually luminous GRBs such as GRB990123. 



3. Dust Destruction by Thermal Sublimation 

As discussed in §3.3 below, in dense regions radiation with /iz^ > 13.6 eV will largely go into 
photoionizing H and H2, and photons in the 7.5-13.6 eV range will primarily be absorbed by H2, 
which is rovibrationally excited by ultraviolet pumping (Draine 1999b). We therefore first discuss, 
in §3.1 and §3.2, sublimation of dust grains under the assumption that dust heating is dominated 
by 1 eV < /iz^ < 7.5 eV photons, and discuss the contribution of 7.5-50 eV photons to dust 
sublimation in §3.3. Using equations (^) and (0), we find that for typical GRB parameters, the 
prompt luminosity in the 1-7.5 eV range is 2 x 10^^ erg/s and the 7.5-50 eV prompt luminosity 
is ~ 5 X 10^^ erg/s. Since our simple analysis overestimates the flux of GRB990123 by a factor 
~ 2, we will take the typical 1-7.5 eV luminosity to be Li-7.5 = 10^^L49erg/s, and the 7.5-50 eV 
luminosity to be Ly.s-so = 2.5 x 10^^L49 erg/s. 



3.1. Grain Heating 

Consider a grain of radius a located at a distance r from a transient source of radiation 
radiating a 1-7.5 eV power Li_7.5(Ar2/47r) into a solid angle If the radiation from the GRB is 
"beamed" into AQ < 47r, we will consider only dust grains within the beam. Note that we expect 
the beaming of optical-UV emission from the forward shock to be similar to that of gamma-ray 
emission. 

The grain temperature T is determined by 

e-^^^Quvvra^ = (Q)T47raVr^ - Attu^^^B , (8) 
47rr^ at m 

where r is the effective optical depth for attenuation of the optical-UV flash, p is the density of 
the grain material, m is the mean atomic mass, B is the chemical binding energy per atom, 

is the usual Planck-averaged absorption efficiency, and Q\jy is the absorption efficiency factor 
averaged over the 1-7.5 eV spectrum of the optical-UV flash. For the grain radii a ^ 10^^ cm 
expected in dense clouds, we will assume Quv ~ 1- Since we are interested in energy depositions 
large enough to sublime grains, the heat capacity of the grain has been neglected in equation (0). 
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The sublimation rate can be approximated by 



(10) 



Guhathakurta & Draine (1989) have estimated i/q ~ 2 x 10^^s"\ B/k = 68100K for Mg2Si04, 
and z/Q = 2 X 10^'^s-\ B/k = 81200K for graphite. We adopt vq = I x 10^^s"\ B/k = 7 x lO^K, 
and p/m ^ 10^^ cm~'^ as representative values for refractory grains. If we assume the grain 
temperature T is approximately constant over the time At = At°^' /{I + z), then the condition for 
the grain to be completely sublimed during this time would he T > Tc, where 



B/k 



In 



2300K 



1 + 0.033 In 



/ a-5 



(11) 



where a_5 = a/10 ^ cm. Equivalently, the grain survival time at temperature T is 



''surv 

(T) 



\da/dt\ 



7.7a_5 exp 



7 X lO'' K 



1 



1 



T 2300 K 



(12) 



The infrared emissivity is quite different for graphite and silicate materials (Draine & Lee 1984; 
Draine 1999a). For the temperature range of interest for dust sublimation, 2000 K ^ T ^ 3000 K, 
we approximate 

^a_5(r/2300K) 



{Qh 



(13) 



1 + ^a_5(r/2300K) 

with A 0.03 and 0.3 for astronomical silicate and graphite, respectively. We adopt equation (^ 
with ^ = 0.1 as representative of refractory grain material. 

Radiation and sublimation then contribute equally to the cooling at a temperature rr=s 
determined by 



T 



2820 K<^ 1 - 0. 040 In 



a_5(rr=s/2300K) 



( Tr=s 

V2800Ky 1 + 0.1a_5(ri.=s/2300K) 



-1 



(14) 



For T > Tj-=s, so that sublimation cooling dominates over radiative cooling, the grain 
temperature is 

Tsub ~ 3030K [l + 0.043 ln{Quye-^ L^gr^^^ , (15) 

where r = lO^^rig cm. 

For r < Tj-=s radiative cooling dominates and the grain temperature is 



'■ rad 



(Q)T/0.iy V 4 

where e~'^ is the instantaneous attenuation of the 1-7.5 eV flash by intervening absorption. 



(16) 
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If attenuation by intervening material can be neglected, grains are completely sublimed out to 
a destruction radius = Rc where Rc is the radius where the unattenuated flash can heat grains 
to the critical temperature T^.. For our nominal parameters, Eq. (11) and (^) give < Tr=s^ so 
that radiative cooling dominates at r = Rc, and the critical distance Rc is given by 

f Li_7.5 guv , , ^, ^^19 ^ QuvL49(l + 0.1a-5) ^'^' 

In the optically-thin limit, then, the optical-UV flash from the GRB will destroy dust in the beam 
out to a substantial distance. 



3.2. Effects of Dust Optical Depth 

While the 1-7.5 eV photons from a GRB may be capable of destroying dust out to Rc ~ 10 pc 
in the optically-thin limit, in dusty regions (such as molecular clouds) attenuation of the radiation 
by dust grains before they are destroyed will limit the grain destruction to a smaller radius. The 
dust optical depth is a function r(r, t) of both space and time, as the GRB flash "burns" its way 
through the cloud. 

The attenuation of ~l-7.5 eV photons is dominated by dust. To estimate the effect of 
high optical depth, we will assume that only a single dust type is present. Let R^ be the dust 
destruction radius: all grains at r < i?^ are destroyed by the heating effects of the optical-UV 
flash. 

We approximate the 1-7.5 eV emission from the GRB as a rectangular pulse. At radii r < i?d; 
the leading edge of the optical-UV pulse is attenuated by the dusty medium through which it 
propagates, but the trailing edge of the pulse is unattenuated since it propagates through a 
dustless medium, and we are neglecting the effects of gas-phase absorption. Rather than solve for 
r(r, t), we will simplify the problem by assuming that the effects of extinction can be approximated 
as primarily a narrowing of the optical pulse, retaining a rectangular proflle. The problem then 
reduces to determination of a function f{r), the fraction of the flash energy which is absorbed by 
dust interior to radius r, and survival of grains at radius r when irradiated by a radiation field 
-^i-7.5/4vrr^ for a time (1 - f )At. 

The function /(r) then satisfies 

df ^ 2 min[tsurv,(l-/)At] 

^ = QvvridTTa — (18) 

where tsurv('") is the survival time of a grain irradiated by the unattenuated radiation field. 
In the unattenuated portion of the optical pulse, grains are heated to temperatures given by 
equations (|l5,16) with r = 0. The dust destruction radius i?d is then determined by the condition 



tsurv(r(i?d)) = [1 - f{Rd)]At. 
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Figure |l| shows Rd as a function of cloud density nn = n(R) + 2n(H2) for a standard 
dust-to-gas ratio nii{4:-K/3)a'^p/niimii = 0.01, for several different values of the GRB 1-7.5 eV 
luminosity (-L49), duration (Ati), and characteristic dust grain radius (a_5). The radius Rd for 
"typical" GRB parameters is shown as the heavy curve. We see that if a GRB occurred in a dusty 
region, the optical-UV flash from the GRB can clear out a substantial amount of dust which lies 
in the beam. 



3.3. H and H2 Absorption 

Radiation with hv > 13.6 eV may ionize H and H2. Neglecting dust opacity, the prompt flash 
will be able to photoionize a mass 

Mion ~ T^TTt ^ 4 X lO^M0L49Ati (19) 

25 eV 

or out to a radius 

Rion ~ 1.2 X 10^^ cm — — - (20) 



where ns = nn/lO^ cm~'^, and we approximate the flash by a rectangular pulse of duration 
At = lOAti s. We have taken the typical 13.6-50 eV luminosity to be ~ 2 x 10^^X49 ergs s~^. 

The dust destruction radius, Rd, is compared with i?ion in Figure |l[ At densities for which 
Rd > Rion, typically nu ^ 10^ cm~^ (characteristic of star-forming regions), our neglect of 
absorption of hi' > 13.6 eV photons by dust in estimating Rion is justified, because the dust will be 
destroyed by absorption of hi^ < 7.5 eV photons prior to arrival of most of the ionizing photons at 
a given point in the cloud. At these densities, hu > 13.6 eV photons are fully absorbed by the gas 
in a region smaller than the dust destruction zone, thus justifying our neglect of ionizing radiation 
when estimating R^- Since the number of 7.5-13.6 eV photons is comparable to the number of 
hu > 13.6 eV photons, this also shows that, for densities nn ^ lO^cm"^, the 7.5-13.6 eV photons 
will be mainly absorbed by H2, and can be neglected for purposes of dust destruction. 

At lower densities, where Rd < Rion, some fraction of hv > 7.5 eV radiation would also 
contribute to dust sublimation. At these densities the grain temperature near Rd is determined 
by radiative cooling and therefore Rd oc L^^'^, where L is the UV luminosity available for dust 
heating. Since Ly.s-so — 2Li_7.5, the contribution of > 7.5 eV radiation may increase Rd by 
up to ~ 50%. 



4. Electrostatic Disruption? 



Because of the large fluence of energetic photons, dust destruction could, in principle, also 
result from extreme ionization of the dust grain. High degrees of ionization could result in fission 
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of the grain, or emission of individual ions from the grain surface by the process known as "ion 
field emission" . 



4.1. Coulomb Explosions? 

For an approximately spherical grain of radius a, charged to a potential U , the tensile stress 
averaged over a cross section vra^ is S* = {U jaf' j^m. If the maximum tensile stress which the grain 
material can support is Smax, then the potential gradient and grain charge Z are limited by 

(!)<^-°^x^°'( io.°/;:o.n-0 "'^"""' 

We note that if a grain with S = S'max fissioned into two halves, each of the fragments would have 
S « 2~^/^5max and there would therefore be no further fragmentation unless additional ionization 
took place. 

A grain contains ~ 3 x 10^0^5 electrons. If the mean atomic number is ~ 10 and the 
photoionization cross section per electron is a « 10^^^ cm^, then substantial grain destruction by 
this process would require a fiuence F{hv > lOkeV) ;^ 2.5 x lO-'^^alg cm~^. Grain fission would 
therefore occur within a radius 

i?Fis ~ ^^vcEl^^all (23) 

where we have taken the fiuence F{hv > lOkeV) (£?/20keV)/47ri?^. Note that this process 
changes the grain size distribution and therefore affects the optical extinction curve, but grain 
fission alone would not appreciably reduce the ultraviolet extinction, and might even increase it. 

4.2. Ion field emission 

Ideal materials have tensile strengths S^ax ~ lO^^dyne cm~^, so that a Coulomb explosion 
would not take place until the surface electric field reaches [//a ~ 3 x lO^Vcm"^. However, in the 
laboratory electric fields exceeding ~ 3 x lO^Vcm^^ are observed to result in "ion field emission", 
where individual ions are emitted from the sample (MuUer Sz. Tsong 1969). As a result, if the 
grain tensile strength Smax ^ lO^^dynecm"^, intense irradiation by hu ^ lOkeV photons will first 
cause the grain to charge up to Zmax ~ 2.1 x 10^([//a)inax, and each subsequent ionization will 
result in emission of an ion (assuming that electron capture is negligible during the ~ 10 s of the 
gamma ray burst). 

If the mean atomic number is ~ 10 and the photoionization cross section per electron 
is a ^ 10~^^cm^, then substantial grain destruction by this process would require a fiuence 
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F{hiy > lOkeV) ^ 10^^ cm ^. Grain destruction by ion field emission would therefore occur only 
within a radius 

i?iFE ~ O.Spc^g/' (24) 
Ion field emission is evidently much less important than sublimation. 



5. Dust Infra-Red Emission 

For clouds of mass M < W'^Mq extending to r ^ 8 pc, most of the optical-UV energy 
absorbed by sublimated dust is re-radiated in the infra-red, typically around A = 1/xm. This is 

1/2 

due to the fact that at distances larger than Rr=s ~ 8LJ^ pc radiative cooling of dust grains 
dominates over sublimation cooling [see Eqs. dHJll)], and to the fact that for M < lO'^M© grains 
are heated to the temperature Tc ~ 2300 K required for complete sublimation out to a distance 
Rd 10 pc [see Eq. (11), Fig. 1] (For M > W'^Mq, optical photons are completely absorbed at 



distances ^ 10 pc, where sublimation cooling dominates and only a small fraction of the absorbed 
energy is re-radiated). 

At radii where radiation cooling dominates, the grain temperature drops approximately as 



T oc r~^/^. The strong dependence of grain survival time tgurv on temperature, Eq. (12), then 
implies a sharp transition, i.e. over a distance Ar <C Rd, between the region at r < R^ where tsurv 
is much smaller than the flash duration At, to the region at r > R^, where tsurv ^ At. Since 
the energy radiated by grains is proportional to max(tsurv; At), infra-red emission of sublimated 
dust is dominated by emission from grains just outside Rd- If the optical depth for UV photons 
due to dust at r > i?d ~ lOpc is ruv ^ 1; then most of the flash energy would be absorbed by 
dust and re-radiated in the infra-red. Most of the infra-red radiation would escape the cloud, and 
may therefore be detected, if the infra-red optical depth is not high, tjr ^ 1. For Qahs,i^ oc v^, the 
requirements ruv ^ 1 and tjr ^ 1 may be written as 0.2 ^ ry ^ 2. 

In order to estimate the dust infra-red luminosity in the case where tir ^ 1 and ruv ^ 1; let 
us first assume that optical-UV flash emission is beamed into a small solid angle around the line 
of sight, Ar2 = ttO'^ <C 47r. The observed duration of the infra-red emission is then AtiR ^ Re^/2c, 
where i? ~ 10 pc is the radius out to which grains are heated to ~ 2300 K. This may be written 
in the form AtiR = 2 (i?/c)( AQ/47r) 20(Ail/0.01) day, which is valid for An = 4tt as well as for 
AQ <C 47r. The infra-red luminosity is given by the ratio of the flash energy absorbed by dust, 
^1-7.5 ~ 10^°(AJ7/47r) erg, and the observed duration AtjR, Ljr Ei^-j,^/ Atm ~ W^^erg/s. 



6. Implications 



The luminosity of the prompt optical-UV emission accompanying GRB 7-ray emission is 
given by Eqs. (P) and (^). For typical GRB parameters we expect an optical-UV flash with 
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1-7.5 eV luminosity Li_7.5 ~ 1 x lO^^erg/s, assuming isotropic emission. Such a UV flash can 
destroy dust by sublimation out to an appreciable distance, i?d ~ 10 pc (see Figure |l]), and may 
clear the dust out of ~ lO''(Ar2/47r)M0 of molecular cloud material, where AQ is the solid angle 
into which the optical-UV emission is beamed, and where dust is sublimed. If GRB sources indeed 
lie in dusty regions, then the extinction would decrease with time during prompt optical-UV 
emission, over tens of seconds. Detection of such time dependent extinction would provide strong 
constraints on the GRB environment. The destruction of dust implies that existing, or future, 
observations of not-heavily-reddened fireballs are not inconsistent with GRBs being associated 
with star formation. 

We have shown in §5 that if the optical depth due to dust beyond ~ 8 pc is of order unity, 
most of the UV flash energy is absorbed and re-radiated in the infra-red, typically at A ~ 1/xm. 
The resulting infra-red luminosity, Ljr ~ lO^^erg/s, extends over an apparent duration of 
~ 20(1 -|- z)(Ail/0.01) day. For GRBs at z = 1, therefore, K-band photometry may reveal thermal 
emission from dust grains. 



In fact, such emission may already have been observed in GRB970228 (Fruchter et al. 1999) 



and GRB980326 (Bloom et al. 1999). In both cases, a deviation from a power law decline of 
optical flux, which at early time is consistent with synchrotron emission from shock accelerated 
electrons, is observed at ~ 30 d delay. As the flux drops below ~ l^Jy, a new infra-red emission 
component is revealed, with a flux fi, ~ 0.5/iJy between A = 2^m and A = l//m for GRB970228 
and fu 0.7^Jy at A = 0.9/um for GRB980326. In both cases, the spectrum is modified at this 
time to fp oc A^ at 0.5//m ^ A ^ 0.9^m. The infra-red flux is of the same order of magnitude 
estimated for dust grain emission, and the spectrum is consistent with dust emission peaking at 
~ Ifim (at the source redshift), provided GRB980326 is at redshift z ~ 0.4. We note that for the 
typical parameters adopted in this paper, dust emission is expected to peak at somewhat longer 
wavelength, ~ 1.5/im. However, since the grain properties are not well known, dust emission can 
not be ruled out as an alternative to the proposal that the "excess" emission is due to a supernova 



(Bloom et al. 1999, Reichart 1999, Galama et al. 1999) We note also that the non-detection 



of optical emission from GRB980326 at t ^ 200 d, implies, under the dust emission hypothesis, 
beaming of the optical-UV flash to 0.01 ^ AQ ^0.1, consistent with the interpretation that the 
rapid, flux decline is due to the flreball being a jet of small opening ang le ( iRhoads 1999| ). 

We have shown in §2 that strong optical-UV emission requires, like GRB 7-ray production, 
large initial expansion Lorentz factor, Fj > 100, which also implies that the plasma emitting the 
optical-UV flash must be expanding with F > 100. Thus, if the flreball is a jet of finite opening 
angle, 1 ^ 9j > 1/F ~ 0.01, then both 7-ray and optical-UV emission will be confined to a 
small solid angle AQ = tt^J. In this case, dust would be evaporated only within a narrow cone 
around the line of sight. A jet-like fireball expands as a conical section of a spherical fireball, with 
F oc as long as F > 1/6 j. After deceleration to F < l/Oj, the jet expands sideways, its 

opening angle increasing to ~ 1/F and F oc ( [Rhoads 1999| ). At this stage radiation reaching 

us must travel a distance / « rsin(F^-'^ — 9j)/ sm{9j) through gas which was not exposed to the 
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initial flash, and which will still contain dust. For 1 ^ F ^ using t ~ r/2r^c, we have 

/ ~ r/T9j ~ 2Tct/6j ~ 0.1pc(t/l week)^/^/^^''^. Thus, a significant increase in extinction over time 
would be observed if Oj <C 0.1. 



Confinement of 7-ray and strong optical- UV emission to a small solid angle is not limited 
to the case of a jet-like fireball. It may also arise if Fj, the initial expansion Lorentz factor, 
is anisotropic. Consider a fireball carrying similar energy per unit solid angle in all directions, 
with Fj a decreasing function of angle with respect to the line of sight, such that Fj <C 300 for 
9 > 6-y. In this case, optical-UV (and 7-ray) emission would be suppressed at angles 9 > 9j. The 
isotropic fireball energy per unit solid angle implies that, after a transition phase, the fireball 
would approach spherical expansion, with F oc t~^^^. At this stage, most of the radiation detected 
at time t by a distant observer is produced by fireball plasma within a narrow ring of radius ~ r/F 
around the line of sight, where the fireball radius r and expansion Lorentz factor F(r) are related 
to t through t = r/2F-^c ( Waxman 1997c| ). Thus, here too a significant increase in reddening is 



expected once F drops below 1/0^. For 1 ^ F <C 1/^^ most of the radiation reaching us must 
pass through a path length / ~ r/F^^ ~ 2Tct/9.y ~ 0.05pc(t/l week)^/^/^^ of gas which was not 
exposed to the initial flash. 

Finally, it should be noted that although optical flashes are expected for typical GRBs [with 
S53 ~ 1 and riQ ~ 1, see Eq. (^)], prompt optical emission has been observed for only 1 GRB to 
date and therefore may not accompany all GRBs. Indeed, fireballs in low density environments 
with n <C 1 cm~^ would not be expected to produce strong prompt emission [see Eq. (P)]. In 
addition, if the fireball initial Lorentz factor Fj ^> 300, reverse shock emission may be shifted to 
photon energies above 7.5 eV [see Eq. where most photons are absorbed by H and II2. 
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Fig. 1. — Radius Rd out to which grains arc destroyed by thermal subhmation, as a function of 
cloud density nn, for different values of L49, Ati, and a_5. The heavy curve is for "typical" GRB 
parameters. The lines Ay = 10 and Ay = 30 show the radius of a cloud having Ay = 10 and 30 
from center to edge. Also shown (broken line) is the radius Rion out to which the gas is photoionized 
by a flash with L4gAti = 1. Dotted lines indicate radii with enclosed gas mass from 10^ Mq to 
1O«M0. 



